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The effect of filtering on the Signal-to-Noise Ratio (SNR) of a coherently demodulated 
band-limited signal is determined in the presence of worst-case amplitude ripple. The 
problem is formulated mathematically as an optimization problem in the L 2 -Hilbert 
space. The form of the worst-case amplitude ripple is specified , and the degradation in the 
SNR is derived in a closed form expression. It is shown that when the maximum passband 
amplitude ripple is 2 A (peak- to -peak), the SNR is degraded by at most (1 - A 2 ), even 
when the ripple is unknown or uncompensated. For example , an SNR loss of less than 
0. 01 dB due to amplitude ripple can be assured by keeping the amplitude ripple to under 
0.42 dB. 


I. Introduction 

Amplitude ripple is inherent in most physically realizable 
filters. We seek to determine the system performance degrada- 
tion resulting from amplitude ripple. The Signal-To-Noise 
Ratio (SNR) is used as the performance criterion. Our results 
provide an easy method to determine the worst-case loss due 
to amplitude ripple. To derive the expression for the worst- 
case loss in the SNR, the class of worst-case amplitude ripple 
is explicitly found. 

We consider filtering a signal y(t) which is composed of a 
bandlimited signal s(t ), with bandwidth (/ 0 — W, f 0 + W) 
Hertz, added to a noise process (rt(f)}, where f 0 denotes the 
center frequency and W is the half bandwidth of the signal. 
The passband of the filter covers the same band of frequency 
as the signal s(f). This filter is shown in Fig. 1. Ideally, the 
transfer function for the filter is 

X 0 (joo) = e^ r 


for Icol e (2rr(/ 0 - IV), 27 t (/ 0 + IV)) where r is a constant group 
delay. The shape of the filter response at frequencies outside 
of the signal band is shown to be immaterial. 

The amplitude spectrum of the ideal filter, H Q (jto), has con- 
stant gain and group delay in the passband frequency. This is 
never really true in practice where the gain of the filter ex- 
hibits a bounded ripple in the passband. The amplitude ripple 
is the deviation of |// 0 (/w)| from ideal and is denoted by 
A(/co). Characteristics of both an ideal and a nonideal band- 
pass filter are depicted in Fig. 1. The transfer function of the 
nonideal filter in the passband is represented as 

= (1 + A(/W» f (1) 

for 1 gl> I e (27 t (/ 0 - W) , 27r(/ 0 + W)) where 0(/co) represents 
the deviation in phase from constant group delay, r. 

A critical issue in the specification of the filter h(t), for the 
design engineer, is to determine the impact of A(/co), and 
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<p(ju>) on the SNR. The degradation of SNR for phase devia- 
tion 0 (/cj) has been studied previously. J. Jones (Ref. 1) in 
1972 analyzed the filter distortion effects of the phase non- 
linearity for BPSK and QPSK, and has shown that when phase 
deviation <H/co) is bounded by </> max in absolute value, the SNR 
is degraded by at most a factor of cos 2 0 max . 

In this article, it is shown that if the amplitude ripple is 
bounded by A (i.e., IA(/cj)K A), the SNR in the presence of 
the worst case amplitude ripple waveform is degraded by at 
most (1 - A 2 ). This result holds even when the amplitude 
ripple is unknown, or known but not compensated. 

To summarize the outline of the rest of this article, Sec- 
tion II describes the system under study and the underlying 
assumptions for which the SNR figure is analyzed. In Sec- 
tion III, the closed form SNR expression is derived for a 
coherently demodulated signal which is filtered by the non- 
ideal transfer function characteristic defined in Eq. (1). In 
Section IV, the class and properties of the worst-case ampli- 
tude ripple are specified. Finally, in Section V, we make some 
concluding remarks based on our results. 

II. Formulation 

We consider a received waveform containing signal and 
noise, that is ,y{t) = s(t) + n(t). The signal amplitude spectrum 
S(jco) is band-limited to I/I e [/ 0 + W, / 0 - W ] , and its wave- 
form is completely known during each te [0,7"] . The noise 
process (rt(f)} i n our analysis is assumed to be an Additive 
White Gaussian Noise (AWGN) process with single-sided spec- 
tral density N 0 W/Hz. The results generalize to the case where 
the noise is not white. The only restriction is that («(/)}be a 
wide-sense stationary process. This implies that it has zero 
mean and autocorrelation function R n (j) - E[n(t) n(t + r)] , 
where E [ •] denotes the expectation operator. 

The optimal receiver for the observed signal y(t ), which 
maximizes the SNR, is a matched filter (Ref. 2). This solution 
is expressed in the form of the Fredholm integral equation of 
the first kind. There are known methods to solve this integral 
equation explicitly to find the optimal matched filter solution 

wo. 

For an AWGN channel the matched filter solution is 
h MA (t) = s(T - t ) or, equivalently, in the frequency domain it 
is H MA {jcE) = S*(/co) q-J^ t . (Throughout the article, super- 
script * denotes complex conjugate while a midline * denotes 
convolution). 

For the case in which the noise is only wide sense station- 
ary (not necessarily AWGN) with spectrum S n (ju>), the 
matched filter solution may be expressed under certain 


assumptions (Ref. 2) as H MA (/a>) - S*(jco ) e^ u>T /S n (/co). 
This transfer function is recognized as the matched filter 
transfer function for the white noise case divided by the actual 
power density of the noise. Therefore, it is possible to gen- 
eralize our result for wide sense stationary noise processes by 
simply using a matched filter which is matched to both the 
noise and the known signal. Thus, with no loss of generality, 
in our subsequent analysis we assume that the noise is white 
and Gaussian. 

In digital communication systems, the signal s(r) is modu- 
lated at the transmitter to a Radio Frequency (RF) by multi- 
plying s(t) by the carrier signal cos(co 0 f), where 2; r co 0 is the 
carrier frequency. At the receiver (Fig. 2), the observed signal 
is filtered and then demodulated by multiplying the observed 
signal by 2 cos(co 0 r). This signal is passed through a zonal 
low-pass filter to filter out the double frequency terms pro- 
duced by the multiplication operation. The output of the 
low-pass filter is then fed to the matched filter. 

The sampled output of the matched filter each T s is de- 
noted by M ( . 

The SNR is defined as the ratio of the square of the ex- 
pected value to the variance of the random variable In the 
following section, a closed form expression for the SNR is 
derived. This expression is formulated in the form of a func- 
tional. We minimize this functional over the ensemble of all 
possible amplitude waveform ripples in the passband of the 
filter h BF (t), as shown in Fig. 1. 

III. Signal-to-Noise Ratio Expression 

Since the filtering processes are linear, we can consider the 
system response to signal and noise separately. We need to 
determine the mean value of the signal and the variance of the 
noise, both at the matched filter output. We denote the 
response of each stage of the system to the signal by e.(r), 
as shown in Fig. 2. The signal is represented by amplitude 
spectrum throughout the following analysis. Throughout this 
article, the square of a complex function is meant to be the 
magnitude square of that function. 

Neglecting the noise response, the amplitude spectrum of 
jy(f), T(/co), is expressed as 

Y(jcj) = S(j(o j - w 0 )) + SO (w + w 0 )) (2) 

We denote the bandpass filter H\jui) as 

H'(joj) = H bf (j(g> - w 0 )) + + w 0 )) 
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where H BF { •) is the transfer function of a complex low-pass 
filter. The spectrum of the output of the bandpass 

filter under study (Fig. 2), is 

E^fco) = //'(/co) Y(jcd) (3) 

The demodulated waveform e 2 (t) has the spectrum containing 
the sum and difference frequencies, 

E 2 (/gj) = (/(w - w 0 )) + E x (/(a> + W 0 )) 

= //'(/(« - co Q )) y(/(w - w 0 )) 

+ #'(/(« + w Q )) r(/(w + co 0 » (4) 

The demodulator is followed by an ideal low-pass filter, which 
filters out the double frequency terms, and the resulting out- 
put spectrum £ 3 (/co) is E 3 (jco) = E 2 (ju)) // LP (/co), where 

( 1, 1(0 I <27 T(f 0 -W) 

H lp (/<•» = 

( 0, otherwise 

Thus, the output of the low -pass filter from Eq. (4) is 

E,(M = (5) 

Note that if the low-pass filter is not ideal, its deviation from 
ideal should be included in the filter under study. 

Neglecting deviation in H MA ( ;co) from ideal, the matched 
filter that maximizes the SNR has the transfer function 
(Ref. 2) h MA (r) = s(T - r), or equivalently 

H MA (ju) = S*(jco)e-^ T ( 6 ) 

Let for simplicity H(j co) denote //^ F (-/co) + H BF (j<jS). Then 
the matched filter output e 4 (f) can be represented as 


e 4 (T) = - [ - jH(fu) \S(ju)\ 2 dw 
1 

The system noise response is denoted by z(f). The random 
process {M ( } (taking values in ir 1 ) is the sum of the filtered 
signal plus the filtered noise component. Hence, we can write 

M. = e 4 (T) + z(D (8) 

Taking expectation of Eq. (8), and noting that the noise is 
assumed to be zero mean, we get 

E[M.] = JL Jh(/co) |S(/w)| 2 dev (9) 

To compute Var [M ( ] , let z(r) = n(t) * x(r). Note that z(t) is 
a filtered white noise process, which is filtered by the filter 
under study and the matched filter. The cascaded filter is 
denoted by x(t). 

Let X{j(jS) = S(/gj)//(/co) e~ /u,7 \ thus we have 

Var [AT] = Var [z(t)\ I (10) 

1 I t=T 

and 

z(r) = 

(ID 

E[n(t)n(t + r)] = N 2 6(r) 

From Eq. (1 1), the variance of z(t ) can be expressed as 

N o C 

Var[z(f)] J\X(jev )\ 2 dev (12) 

/ 

Combining Eqs. (12) and (10) and evaluating these expressions 
at t = T, results in 


£ 4 (/w) = £' 3 (/w)5*(/co)e‘ / “ r 


and substituting E 3 (ju>) using Eq. (5), e 4 (r) may be expressed 
as 


e 4 (0 = ^ Jnijev) |S(/w)| 2 dco (7) 

where / = [-H/W] . The output of the matched filter is sampled 
at the end of every time interval T. Thus at t = T we have 


Var [/If.] 


M - * fl 

< J 4tt J 1 


|//(/co)5(/co)| 2 dec 


Therefore since SNR = (/^[AfJ ) 2 / Var[Af.] we have 


snr = - 


^H(jev)S 2 (]ev)dev^ 
o J' |//( /gj) S(Jev) | 2 dw 


(13) 


(14) 
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The filter transfer function is //(/co) = H BF (j<jS) + H* f {-]<jS) Using the result of Theorem 1, we can state the following 
and furthermore H(jco) is corollary. 


H(jio) = (1 +A(/co))e /OJT VM<2 nW 

We assume that the group delay is negligible (i.e., e ^ T 1), or 
it is compensated in the matched filter, hence, by substituting 
H(jco) in the SNR expression (14), the SNR is 


Corollary 1. Define interval /,, as any subset of the inter- 
val / such that 

J S 2 (j<x>)du = 


snr -^ 


J(\ + A(joj))S 2 (/aj)du 2 
J(] + Af/co)) 2 S 2 (jaj)du 


(15) 


Then Eq. (15) is minimized by 

! A, for ue/j 
-A, for o> ^ / 


IV. Worst Case Amplitude Ripple 

With no loss of generality assume / = [0,1 ] and let 

£ 

J S 2 (joj)dco 
o 

Using Eq. (15) we can formulate the minimization problem 


Corollary 2. The minimum SNR for Eq. (15) is achieved by 
A(/co) of the form specified by corollary (1), and further- 
more, the total SNR for the worst case ripple distortion is 
(1 - A 2 ) SNR ideaI . 

Corollaries (1) and (2) are direct consequences of Theorem 
1. The point a , in Fig. 3, indicates the point at which the inte- 
grated power of the signal s(f) is (1 - A)/2. 


(I 

I 


(1 +*©)</{ j 2 

F(g) = inf A_J1 Z_ 

[ 0,1 ] 


(1 + g(H)Y d% 


subject to the constraint 

l^(f)|< A< 1 (16a) 

The integrals are understood in the Lebesgue sense. 


In the appendix we prove the following theorem. In the 
proof, the class of the functions for which the minimum 
occurs is explicitly exhibited. 


In general, the SNR is minimized when the ripple is +A for 
frequencies containing (1 - A)/2 of the energy of s(t), and 
-A for frequencies containing (1 + A)/2 of the energy of 

s(t). 

The shape of the amplitude ripple is not unique, and it is 
the whole continuum of step functions which satisfies the con- 
dition stated in corollary (1). To construct another amplitude 
ripple waveform which satisfies the conditions of corollary 
(1), one can take the waveform of Fig. 3 and move a segment 
from [0,tf] to (tf,l] , and move an equal energy segment from 
(<z,l] to [0,ff] . Conceptually, this method may be thought of 
as juggling equal energy line segments from each interval. This 
method results in obtaining a new step function from the 
basic function of Fig. 3 which satisfies the statement of 
corollary (1). 


Theorem 1. There is a continuum of measurable step func- 
tions which minimize Eq. (16), and at the minimum F MIN (g) = 
1 - A 2 , and|X£) is 

! A, for 0 < £ < 

-A, for < | < 1 


V. Discussion and Conclusion 

We have shown that the worst-case loss in SNR due to 
amplitude ripple is (1 - A 2 ). This result can be used to specify 
filters for communication receivers, such as the Advanced 
Receiver for the NASA’s Deep Space Network. 

To express the maximum allowable ripple, for a given loss 
(Z, dB ) in dB, we expand the expression for dB of ripple and 
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keep the first terms. Thus for small ripple and loss, the ripple 
loss in dB may be expressed as 

Ri PPl e dB ( L d B ) = 4168 ' /Z ^ 

Consider a simple example: If the system can tolerate a loss 
of 0.01 dB due to amplitude ripple, the maximum allowable 


ripple in dB would be 0.4168 dB. The loss is much less than 
the ripple and decreases as the square of the ripple. 

From a mathematical standpoint, we solved the minimiza- 
tion problem stated in Eq. (16). And we showed that the opti- 
mal solution lies on the boundary, and it is a continuum of 
measurable step functions in the interval [0,1]. This is expli- 
citly exhibited in the appendix. 
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S(t) + n(t) 



y(t) 


A(ja;) -PASSBAND 



Fig. 1. Filtering 



Fig. 2. Communication model 



Fig. 3. Ripple waveform 
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Appendix 


Let 0 < A< 1 and 


We embed C n in C’ n+2 by 


C A = { /: [0,1] * 1R 1/ measurable 

and | /(z) - 1 1 < A for almost all u [ 


( fl 0> • • • > a n 9 C 1 ’ ■ * * ’ c n) 

^(a 0 ,a 0i . .. i a n _ l ,a n ,a n \ 1 + A, Cj , . . . , c n , 1 - A) 


We equip C A with the induced metric from L l (I) (/= [0,1] ) so 
that C A becomes a complete metric space (we identify func- 
tions equal almost everywhere). Consider the continuous func- 
tional F on C A defined by 


nn = 



In the appendix we investigate the existence and nature of 
minima of F, Let 


and prove the assertion for C^ +2 . Let 


M 

y 



fy (x)dx 


fy ( X ) ^ 


Assume f assumes value c i=- 1 ± A on a nonempty open inter- 
val (a i _ 1 , a t ). Either c >M y or c <M y and assume for definite- 
ness that the former possibility occurs. We may assume c 2 = 
c>M y . In fact it is easy to see that there is y' e C n+2 such that 
My’ = M y , F(fy) = F{f y 0 and c 2 = c where y = (a’ 0 , a [ , . . . , 
a n +2 > c i > • • • » c n+ 2 )* Let 0 < 5 < | a 2 - | and define © e C^ +2 
by 


C = | . . . , a ; c, , . . . , c ) | 

n t v 0 ’ n 1 ’ * n J ' 

0 = < a, < . . . < a = 1, c. e [1 - A, 1 + A] [ 

01 n ’ j L J J 

Then is naturally identified with a compact subset of 

[0,1] n “ 1 X [1 - A, 1 + A] n . To every y e C n we assign the 

function f y e C y which takes value c\ on the open interval 
j, a t ) provided it is nonempty. (Note that f is undefined 
on a finite set, and can be assigned arbitrary values there.) 
The mapping 


7 —F(f y ) 

is a continuous function on C n and hence achieves a minimum 
onC„. 


Lemma. Let y = (a Q , . . . , a n \ c v . . . , c n ) and assume F(f y ) 
is a minimum on C n . Then c ( = 1 ± A whenever a i _ l <a r 

Proof. Let C' n+2 be the compact subset of C n+2 consisting 
°f ( a o ’ • • • * * ,c„ +2 ) such that 

Cj = 1 + A and c n + 2 = 1 - A 


0 = (a 0 , a x + M 2 , • • • ,a„ +2 ; 1 + A, c 2 , . . . ,c„ +1 , 1 - A) 
Then by a simple computation 


F (4) 


^4 + a 
£ + 0 


where 


and 



fy(x)dx 


fy (x) dx 


f 

a Jo 


f y (x) dx + 2 (1 + A - c) 5 


c+y(l + A-c) 


Since c > M 7 , by taking 5 > 0 sufficiently small we can ensure 
ot/(t < A/B and consequently F(f & ) < F(f y ). This proves the 
lemma. 
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Therefore inf F(f) is achieved for some F(f @ ) with © = 
(0, a, 1 ; 1 + A, 1 - A). It is a simple exercise to show that a = 
(1 -A)/2 and F{f & ) - 1 - A 2 . 

Finally we note that for any partition / = £,u£ 2 with 
meas^j ) = (1 - A)/2, meas(E 2 ) = (1 + A)/2 the function 

( 1 + A, if x e £j 

/(*) = } 

( 1 - A, if x eE 2 


In view of the lemma we may assume f y takes only values is also a minimum for F. The argument above also shows that 

1 ± A and hence there is 0* e C 2 such that F(f lk ) = F(f &k ). all minima of F are of this form. 
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